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ABSTRACT 



o 
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Typically, amorphous organic materials contain high density of traps. Traps hinder charge transport and, hence, 
affect various working parameters of organic electronic devices. In this paper we suggest a simple but reliable 
method for the estimation of the concentration of deep traps (traps that keep carriers for a time much longer 
than the typical transport time of the device). The method is based on the measurement of the dependence of 
the total charge, collected at the electrode, on the total initial charge, uniformly generated in the transport layer 
under the action of a light pulse. Advantages and limitations of the method are discussed and an experimental 
(N ' example of the estimation of the density of deep traps in photoconductive organic material poly(2-methoxy-5- 

■ (2'-ethylhexyloxy)-l,4-phenylenevinylene (MEH-PPV) is provided. 
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1. INTRODUCTION 

Most organic materials used in nowaday organic electronic devices are amorphous materials (organic glasses) P 
For that reason traps are ubiquitous in solid organic semiconductors. They could be chemical impurities, which 
are intrinsic to the organic material (in many devices organic polymers are used and polymers are notoriously 
difficult to purify) or accumulated as products of the degradation of the organic material during the device oper- 
1 ation, or they could be structural imperfections of the material; the possible nature of traps is very divers. From 
the point of view of the operation of electronic organic devices the most important and common characteristic of 
O ' traps is their ability to capture charge carriers and keep them for a long time, thus hindering the charge transport 
and leading to the degradation of the performance of the devices (solar cells, light emitting diodes etc.jPS In 
some situations traps significantly modify transport characteristics of organic materials, such as the dependence 
t-H | of the carrier drift mobility \x on the applied electric fieldP^ The most dangerous type of traps is the deep traps 
capable to keep the captured carriers on the time scale much longer than the relevant time scales of the device 
(i.e., infinitely). In some cases traps are favorable: consider, for example, a photorefractive process, where one of 
the major steps is the trapping of carriers of one particular signP In this report we suggest a method to estimate 
the concentration of deep traps in photoconductive organic materials and provide an experimental example of 

■ the estimation of the density of deep traps in the typical organic transport material MEH-PPV. 
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(TJ ■ 2. EXPERIMENTAL AND RESULTS 
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In the experiments, two-photon absorption of the laser beam with the wavelength of 1064 nm in a MEH-PPV 

film has been used for the generation of electron-hole pairs. Weak absorption of the light in the film guaranties 

a uniform spatial distribution of generated charges. Electron-hole pairs were dissociated by the applied electric 

H ■ field F and then drifted charges were collected at the electrodes. 
C3 ■ 

The test specimens consisted of a layer of MEH-PPV sandwiched between the ITO and top aluminum 
electrodes. MEH-PPV from Aldrich was dissolved in toluene. The polymer film was deposited onto ITO/glass 
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substrate by the doctor blade technique in air, and then dried in Ar atmosphere during 10 hours at 80°C. The 
thickness of the films varied between 1.5 /iin and 2.0 /im. The standard Nd:YAG laser (pulse halfwidth of 25 
ps) was used as a light source and the electric signals were recorded by Tektronix TDS3032 oscilloscope. The 
experimental set-up is sketched in Fig. [T^,. 

Typical dependences of the total extracted charge Q on the light pulse energy E are shown in Fig. [TJd. For 
large E, where all deep traps are filled, Q is approximately equal to the total generated charge Qo °c E 2 , hence, 
Q oc E 2 . Indeed, the slope of logQ - log-E plot is very close to 20 For much smaller E another dependence was 
detected 

Q oc E a (1) 

with a ~ 3. Certainly, this dependence reflects the filling of traps with movable charges. Point of intersection 
of the two tangents to the linear regions of the log Q - log E plot may be used as an estimation for the total 
trapped charge Q t and, hence, to the total number of deep traps Mo = Qt/e or trap concentration too = Mq/LS 
(here L is the sample thickness and S is the sample area) . If this assumption is valid, then the intersection gives 
m ~ 5 x 10 13 cm -3 . 

The major problem of this approach is a difficulty to motivate Eq. ([T]) with a — 3 for the relation between Q 
and E in the low energy region. In the next section we consider a simple but realistic model of the trap-controlled 
carrier transport and suggest a more robust procedure for the determination of too. In the rest of the paper we 
mostly use the total number of the collected carriers N = Q/e instead of the total extracted charge Q. 
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Figure 1. (a) Experimental set-up. (b) Extracted charge Q versus light pulse energy E for three applied voltages U : 20V(»), 
40V(A), and 80V(B), correspondingly. Corresponding values of the electrode charge CU are indicated at the right side of 
the plot. Capacitance of the specimen was C = 500 pF. For this particular device the relation mo(cm~ 3 ) = 1.25x 10 23 Qt(C) 
is valid. 

3. MODEL OF TRAP-CONTROLLED CHARGE TRANSPORT 
3.1 No-recombination case 

Figure QJ> indicates that the total extracted charge is less than CU, hence, we can assume that the electric field 
in the sample is uniform in space and equal to the applied electric field F. This means that the carrier velocity v 



*For a comparison, in the case of one-photon excitation, where the light pulse is mostly absorbed in a thin surface 
layer, the dependence of Q(E) is much weaker than Q oc E for high E. 



is a constant in the bulk of the sample. As it was noted in Sec. [3J we may also safely assume a spatially uniform 
initial distribution of the generated carriers. 

Let us begin our consideration with the case when carrier recombination is negligible. Then we can consider 
charges of the opposite signs separately and limit our consideration to the charges of one sign only, with concen- 
tration n(x,t). Concentration of the empty traps is m(x,t). Initial conditions are n(x, 0) = hq, 771(2;, 0) = mo, 
so the total number of generated carriers is Nq — hqSL. Dynamics of the model is governed by equations 
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which take into account the drift of charge carriers with average velocity v and a nonlinear term describing 
trapping of carriers (k is the trapping rate constant). We do not take into account carrier diffusion, a brief 
motivation is provided in the Appendix. 

If the initial distribution of carriers is uniform in space, then drift does not produce a spatial variation of 
n(x,t). The only effect of the drift is moving of the rear front of the carrier distribution with velocity v. Hence, 
for x > vt n(x,t) — n(t) and we can solve Eqs. @ and ([3]) without the drift term, while for x < vt n(x,t) = 
and m(x, t) = mix, x/v). 

Taking into account the conservation law n(x, t) — m(x, t) = tiq — tuq = A, we can write 
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Evolution of the carrier and empty trap distributions is shown in Fig. [3] 
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Figure 2. Relative densities of free carriers n(x,t)/no (solid lines) and empty traps m(x,t)/mo (broken lines) for three 
particular moments of time (shown in the plots); q = 1.5, /3 — 0.3. 



The total number of extracted carriers is 
N = S' 
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This formula is the main result of the paper. General behavior of N(no) is shown in Fig. [21 Taking into account 
that A^o oc E 2 , we see that Eq. ([6]) gives the desired relation between Q and 

^Strictly speaking, the relation Q(E) should include contributions from the carriers of both signs, thus the more proper 
analogue of Eq. © should have two separate terms with different parameters q and /3. In future we will not use this 
more strict relation because the limited accuracy of the experimental data gives no possibility for the reliable extraction 
of individual parameters q and /3 for positive and negative charges. 
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Figure 3. General dependence of N on no in the double linear plot (left) and log-log plot (right) for different values of 
/3: 0.01, 0.03, 0.1, 0.3, 0.5, 1, 3, 10, and 100, from the bottom curve upward, correspondingly. Thick line is the limit 
dependence N/Mo = q — 1. 

Parameter ft — itrapAdrift is the ratio of the characteristic trapping time itrap = (fc^io) -1 and carrier drift 
time idrift = L/v. This helps to understand most important features of Fig. [3] If /3 3> 1, then trapping is not 
effective and N « Nq. In the opposite case of effective trapping /3 -C 1 



In fact, Eq. © is a universal asymptotics of N(Nq) for Nq — >• oo and any finite (3, even for (3 ^> 1, because for 
a very high concentration of carriers all traps should be filled. Indeed, ii q ^> max(l, f3), then 



For (3 ^ 1 this asymptotics develops only for a very high concentration of carriers q ^> (3 (for this reason it is 
not seen in Fig. [3] for j3 ^> 1). 

The case of large (3 (fast transport) is very unfavorable for the determination of Too- Thus, in such a case 
the experiment should be performed for thick samples or a weak electric field (ensuring low v), both conditions 
leading to the decrease of (3. 

For g < 1 the concentration of empty traps is approximately constant, thus the only significant Eq. 
becomes linear. For this reason we must have a linear dependence N — a(f3)No. This is indeed the case 



3.2 Contribution from the carrier recombination 

Let us check, how sensitive are obtained results to the possible contribution from the carrier recombination. In 
this case we have to consider two kinds of movable carriers with concentrations n(x, t) and p(x, t). Dynamics of 
the model is described by equations 



N = 0(f3), n < too, 

N = No-M + o(l), n >m . 
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where k r is a recombination rate constant, and we omit trapping for the simplicity sake. Basic equations are 
almost the same as in the previous caseH but the reaction (recombination) zone is different because both species 
are movable. Again, n(x, t) = n(t) and p(x, t) — p(t) for v n t < x < L — v p t. 

Let us limit our consideration to the case no — pa, so n(t) — p(t) (if the equality does not hold it implicitly 
means some sort of trapping with trapped charges avoiding recombination). Temporal dependence of n(t) is 
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and the total number of carriers extracted at the collecting electrode is 
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Finally, 
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We conclude that the recombination is not important for A^o -C 77, while for iVo S> 77 it produces the dependence 
N (x IniVo. Hence, an experimental observation of the dependence N oc E 2 oc Nq for high E (and, hence, high 
initial number of carriers) is a clear indication that the recombination is negligible for the whole range of E. 
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Figure 4. Simulated data for /3 — 0.04 (a) and /3 = 0.3 (b), correspondingly, obtained by adding small random contributions 
to Eq. (0. We used mo = 1 in some arbitrary units in both plots. Plots show only a part of the full curve (for sufficiently 
high initial number of carriers), providing a more close resemblance with Fig. [Tb. Straight lines indicate the power law 
fit TV oc n^ 2 oc E a , and broken lines show the result of the fit of the data to Eq. (JSJ). Parameters f3^ and 
by fitting, agree well with the true parameters. 



'In fact, recombination kinetics may be considered as a special case of the trapping process for q = 1. For this case 
Eq. (0 directly leads to Eq. (|15|l with k r = k and v p — 0. 



3.3 Determination of a trap concentration from the experimental data 

Our consideration shows that for Nq 3> Mq there is a natural relation 

N w N - M cx E 2 - const, (16) 

yet for smaller Nq there is a transition to another linear dependence 

N « a(/3)N a oc E 2 (17) 

(see Fig. Hence, in this model there is no dependence of the kind N oc E a with a > 2. In fact, such 
dependence hardly exists in any model of the trap-controlled transport because the linear dependence of N on 
Nq for small concentration of carriers no <C mo is a very general property of the charge transport, it just follows 
from the approximate constant concentration of empty traps in such a case. Nevertheless, the transition region 
from Eq. (IT71) to Eq. (fH)]) can simulate this very kind of the dependence in a limit range of E (see Fig. [31 the 
right plot). 

If P -C 1, then the concentration of traps could be estimated by plotting of the experimental dependence N 
vs E 2 and drawing the linear asymptotics for the region of large no . If j3 ~ 1 , then the most reliable way of the 
estimation of Too is the direct fit of experimental data to Eq. 

Figure [4] shows the result of fitting of a primitive simulation of the experimental data obtained as 

N sim = N(n a )(l + 5), (18) 

where AT (no) is calculated using Eq. (JBJ) and i5 is a random Gaussian number with zero mean and magnitude 
a = 0.1, which models a noise in the experimental data. One can see that the fit to the power law dependence (JT|) 
produces strongly overestimated value of mo, and the difference between m rue and m^* progressively increases 
with the increase of /3. Fit of the simulated data to Eq. © produces a reasonable agreement. 

Figure [3] indicates that for f3 > 1 the reliable fit of the experimental data to Eq. ^ in the log-log plot is 
hardly possible. Reasonable question is: is it possible to estimate mo and /? in such a case by fitting the data to 
Eq. © in double linear coordinates? Figure [5^ shows that this is possible for mo if experimental errors are not 
too large, but the estimation for j3 is still not reliable. 

At last, we tried to fit the experimental data, presented in Fig. [lb, to Eq. © (see Fig. [5b). Again, true 
concentration of traps is much lower than the one, estimated by the simple power law fit according to Eq. ([1}. 
It is worth to compare Figs. 0Jd and [3b. In both cases for the low E region a ~ 3, and the corresponding ratios 

power law / right fit , , t 

m o / m o are P re tty close. 



4. DISCUSSION: WHEN THE METHOD CAN BE APPLIED? 

Let us consider the applicability of the suggested method to real organic amorphous semiconductors. The 
necessary conditions for the direct applications of the method are: 1) spatially uniform generation of initial 
electron-hole pairs; 2) carrier recombination can be neglected; 3) the total trapped charge Qt — eMo must be 
less that the extracted charge Q — eN for the high end of the spanned range of E; 4) electric field is a constant 
in the bulk of the sample, so the total charge in the sample must be much less than the electrode charge CU, 
and 5) charge transport should not be very fast (f3 < 1). 

First condition means in our case that the absorption of the light in the sample must be weak. The two- 
photon absorption as a way for the carrier generation is not important by itself for the application of the method. 
The only reason to use this particular method of the generation of initial carriers is a weak absorption of the 
1064 nm light in MEH-PPV, thus ensuring the spatially uniform initial distribution of carriers. Other methods 
of carrier generation, capable to produce the uniform distribution, can be used as well. 

In the second condition we mean a non-geminate recombination of charge carriers, because the geminate 
recombination could be phenomenologically taken into account simply by the redefinition of the proportionality 
coefficient c in the relation between the light pulse energy and total number of initially generated free carriers 
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Figure 5. (a) Simulated data for /3 — 3. Broken line is the result of the fit of the data to Eq. JB). (b) Fit of experimental 
data for MEH-PPV, shown in Fig. [TJj, to Eq. ((6]) (solid lines). Concentrations of deep traps, obtained by fitting, are 
shown in the plot, for different values of the applied voltage. These concentrations do not differ significantly, as it should 
be for the intrinsic parameter of the material. Some difference in values of mo for different values of U is not surprising 
taking into account the poor accuracy of the experimental data. 



A*o = cE 2 . Absence of the non-geminate recombination can be reliably established by checking the linearity of 
the relation Q oc E 2 for high E. 

The third condition is a principal one, because if it is not valid, then we are dealing with the linear regime 
of low trap filling described by Eq. (flOl) , where the only relevant trap-related parameter is km® (see Eq. ([2j for 
the case m(x,t) w mo = const). Naturally, the separate extraction of k and mo is not possible in this regime. 
At the same time, this condition is not an obstacle per se for the use of the method, it just dictates the proper 
range of E, most favorable for the determination of the trap density. 

Last two conditions provide more serious restrictions because in some situations they are mutually contradic- 
tive. Indeed, the applied voltage U should be high enough to provide the uniformity of the electric field in the 
sample, but at the same time the carrier velocity v grows with the voltage, thus invariably shifting the sample to 
the unfavorable regime of high (3 oc v. Yet the example of the MEH-PPV device indicates that for very typical 
organic materials the suggested method could be used successfully. 

In the basic Eq. ([2]) we omitted the diffusive term. Short analysis in the Appendix shows that the diffusive 
contribution to N may be safely neglected in most cases. 

5. CONCLUSION 

In this paper we suggested the simple method for the estimation of the density of deep traps capable to keep 
charge carriers in amorphous semiconductors for a long time. One of the most important advantages of the 
method is a simplicity of the experimental set-up, which includes only standard widely used equipment. In the 
previous Section the major conditions, necessary for a successful application of the method, have been discussed. 
It turns out that these conditions are not very difficult to fulfil, thus the method may be considered as an 
almost universal one for the estimation of the density of deep traps in amorphous organic materials. Successful 
application of the method to the estimation of the concentration of deep traps in typical organic semiconductor 
MEH-PPV supports this conclusion. 



APPENDIX. CONTRIBUTION FROM THE CARRIER DIFFUSION 



There is another limitation of the suggested method, directly related to the structure of our transport model: 
there is no diffusion term in Eq. ©, we set the diffusion coefficient D = 0. Let us estimate when the diffusive 
contribution is negligible. Assuming Lv/D ^> f (weak diffusion), the diffusive contribution to N could be 
estimated in the following way. At t = the distribution of movable carriers is uniform in space. It means 
that the diffusion provides the most important contribution at the rear front of n{x, t), the only place where the 
gradient of n is large. The width of the rear front increases with time as A r (t) ~ (2Di) 1,/2 , or, equivalently, 
with the traveled distance x = vt as A r (x) ~ fiDx/v) 1 / 2 . At the location of the rear front the concentrations 
of carriers and empty traps are, correspondingly 

nM = ^rrP 1 — wt ( 19 ) 

uq — mo exp |— kxyriQ — mo)/ v\ 
m r (x) = n r (x) — no + mo. (20) 

Diffusive addition SN to the total N can be estimated as a number of carriers, trapped in the diffusive zone with 
width A r , moving with the rear front of the carrier distribution 

kS f L 

SN ~ — / dxA r {x)n r {x)m r {x). (21) 



We neglect here a variation of n(x,t) in the vicinity of the rear front, this changes the estimation (|2"Tj) by the 
factor ~ 0(1). Parameter A r varies slowly with x, it can be estimated as A r (L) and taken out of the integral. 
Next, we consider only two limiting cases no ^ wo and no <C mo- If no 3> mo, then 

n r (x) ~ n , 

m r {x) ~ mo exp (— kxno/v) , 



SN mo ( L> x 1/2 



N no \Lv 

Estimation (|2"21 is indeed small. For the opposite case mo no 

n r (x) ~ no exp (— kxmo/v) , 
m r (x) ~ mo, 

SN I / D \ 1/2 



[I - exp (-kLno/v)] . (22) 



N 13 \Lv 



(23) 



To obtain Eq. (|2U)) we used Eq. (|TU| . which relates N with Nq for nio 3> "-o- For (3 ~ 1 this is a small 
contribution, but it can be large if /3 <gC 1. Hence, in some cases of a slow charge transport the ratio SN/N may 
be not negligible, but this is most probable for a very far left end of the spanned range of E, because the factor 

1/2 

(D/Lv) ' is typically small in amorphous organic materials at the room temperature. 

This factor may be estimated using the Einstein relation u = eD/kT. This relation is not strictly valid in 
organic amorphous materials}^! but could be used as a very crude estimation of D. According to the Einstein 
relation D/Lv ~ kT/eU . This factor is indeed very small for the typical applied voltage U ~ 10 -j- 100 V and 
the room temperature with kT — 0.026 eV. We conclude that the typical diffusive contribution to Eq. (J6j> is not 
very important. 
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